Testing the broad applicability of the PBEint GGA functional and its one-parameter 

hybrid form 
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We review the performance of the PBEint GGA functional (Phys. Rev. B 2010, 82, 113104) re- 
cently proposed to improve the description of hybrid interfaces, and we introduce its one-parameter 
hybrid form (hPBEint). We consider different well established benchmarks for energetic and struc- 
tural properties of molecular and solid-state systems as well as model systems and newly developed 
benchmark sets for dipole moments and metal- molecule interactions. We find that PBEint and 
hPBEint (with 16.67% Hartree-Fock exchange) yield the overall best performance, working well 
for most of the considered properties and systems and showing a balanced behavior for different 
problems. In particular, due to their well-balanced accuracy, they perform well for the description 
of hybrid metal-molecule interfaces. 

PACS numbers: 71.10.Ca,71.15.Mb,71.45.Gm 



INTRODUCTION 

Density functional theory [l[ (DFT) is nowadays one 
of the most popular computational methods in quantum 
chemistry and solid-state physics. Despite its success 
however, DFT still suffers some limitations due the ap- 
proximations used in the exchange-correlation (XC) func- 
tional, which describes all the quantum electron-electron 
interactions. In fact, the development and optimization 
of advanced XC functionals is currently a very active field 
of research Q . 

In the years numerous different XC functionals have 
been developed, ranging from local density approxima- 
tions (LDA) to the most advanced orbital-dependent 
XC functionals including exact exchange and se- 

lected exact-correlation contributions. The largest pop- 
ularity in practical applications is however own by two 
classes of functionals: generalized gradient approxima- 
tions (GGA) and hybrid functionals. The former are in 
fact, the method of choice for the investigation of large 
systems (e.g. in biochemistry or solid-state physics) due 
to their very favorable cost-to-accuracy ratio, while the 
latter are the workhorse for computational chemistry. 
Moreover, the GGA functionals attract theoretical in- 
terest because they are the building blocks of the more 
advanced meta-GGA, hybrid, and hyper-GGA function- 
als. 

Due to their simple form GGA functionals cannot ful- 
fill all the exact constraints of the XC energy and cannot 
be accurate for both atoms and solids ^ [l^l • There- 
fore, some selective criteria must be employed for their 
construction. Popular approaches are to develop rather 
specialized functionals by fitting to specific problems and 
training sets or by requiring the satisfaction of selected 
exact constraints of the XC energy. As a result, GGA 
functionals for molecules (PBE APBE revPBE 



[il,RPBEjl3, BLYP fTI[3, OLYP [ll,[l3) or solids 
(PBEsol flo|rARPA-F m, AMOS [H, Wu-Cohen 0) 
are obtained. On the other hand, hybrid functionals 
rely strongly on the underlying GGA approximations and 
share most of their fundamental drawbacks. In addition, 
they might suffer from a certain degree of empiricism 
used to construct the partial inclusion of the Hartree- 
Fock exchange. Thus, hybrid functionals generally im- 
prove over the GGA description but cannot provide a 
really homogeneous level of accuracy for a wide range of 
applications, showing occasionally important failures for 
specific problems (e.g. transition-metal chemistry 21| or 
metal clusters pjj). 

Recently, a growing effort was devoted towards the de- 
velopment of GGA functionals that could provide a well 
balanced description of a large number of properties. In 
fact, a wise selection of important exact constrains of the 
XC energy (as e.g. in PBE or APBE) can be used to op- 
timize the performance of GGA functionals in different 
contexts [23| . Two noteworthy examples are the PBEint 



2J] and the HTBS [25| functionals. These functionals 



are not only equally accurate for molecular and solid- 
state energetic and structural properties, but potentially 
important for interface physics, surface science and clus- 
ter chemistry. Furthermore, they constitute an optimal 
starting point for the construction of hybrid functionals 
of broad applicability. 

In particular, the PBEint functional has been devel- 
oped to respect different exact constraints of the XC en- 
ergy and to connect the slowly-varying density regime, 
where the second-order gradient expansion (GE2) of the 
exchange energy is correct, and rapidly- varying density 
regime, where the PBE behavior is accurate. Thus, it 
is reasonably accurate for an important set of molecular 
and solid-solid state properties and, thanks to its well bal- 
anced behavior, especially suited for interface [l^l, metal- 
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cluster '22'] and surface problems. We recall also that 
the PBEint XC hole density shows exact properties [26| . 
beyond the PBE and PBEsol ones 12711 . related to the 
accurate analysis of metallic surfaces [2g . 

In this paper we propose a short review and a more de- 
tailed assessment of the performance of PBEint for dif- 
ferent problems in chemistry and physics. In addition, 
we use the PBEint functional to build a global hybrid 
functional and investigate the results of this for differ- 
ent amounts of the Hartree-Fock (HE) exchange mixing. 
We find that a rather small amount of HE mixing (only 
16.67%) gives the best results over a broad range of prop- 
erties. We recall that global PBE and PBEsol hybrids 
are instead commonly constructed using 25% and 60% 
of HE mixing jl^, respectively. Thus, the PBEint ex- 
change construction is more compatible with the exact 
exchange, a result that was already proved in the case of 
jellium surfaces [26 1. 



THEORY AND COMPUTATIONAL DETAILS 

The PBEint functional is constructed considering an 
exchange functional 



E^^J p{v)er'ip{r))F4s{r))dv, 



with a PBE-like enhancement factor [11 1 



J^x(s) = 1 + K 



1 + /iS^/ K 



(1) 



(2) 



where 



P 



the electron density, e 



unif 



— (3/4)(3/7r)^/'^p^/'^ is the exchange energy per particle of 
the uniform electron gas, and s = | V/3|/(2(37r^)^/^p'*/'^) 
is the reduced gradient for the exchange. In Eq. ([2]) the 
value of the parameter k can be easily fixed to 0.804 
by imposing that the Lieb- Oxford bound [s^l for the 
exchange energy holds locally jll)] (i.e. for the exchange 
energy density at each point in space). The value of the 
fi parameter can be instead fixed by noting that Eq. ([T]) 
has, through second-order, the gradient expansion 



GE2) 



(3) 



so that /i can be just identified with the second-order 
coeHicient of the exchange gradient expansion. However, 
no unique value can be fixed for fj, with this condition. 
In fact, in the slowly- varying density limit (s — > 0) /i = 
nGE2 _ 10/81 (sij, while in the rapidly- varying density 
limit (s > 1) ^ = /i^^^ = 0.2195 was found to be a 
good choice. Finally, from the semiclassical-atom theory 
it was found = h'^GEi ^ o.26 [H [Si. For the PBEint 
functional we therefore abandon the idea of a constant fi 
and consider instead an s-dependent /i with the following 



ansatz 



Ms) = M^^2 + (A."---A.--) 



PBE ,,GE2\ 



as 



1 + as2 ' 



(4) 



with a = 0.197. Equation ^ assures that the following 
minimal constraints are satisfied: i) in the slowly- varying 
density limit /i — > i^^^] ii) in the rapidly- varying density 
limit /i — !■ /i^^^; (iii) the fourth-order term (oc s**) in the 
exchange gradient expansion vanishes. 

Due to Eq. (|3|) the exchange enhancement factor of 
PBEint varies smoothly between that of PBEsol (at small 
s values) and that of PBE (at larger s values), granting 
a good description of all possible density regimes and of 
many different systems from molecules to solids. Note 
however that the PBEint functional (unlike HTBS) is 
not a simple interpolation between PBE and PBEsol, but 
provides instead a physically meaningful (although inter- 
polated) value for the /i parameter at different density- 
regimes. In fact, a similar construction proved to be 
useful also in_the case of noninteracting kinetic energy 
functionals ^33|. 

For the correlation part, the PBEint functional utilizes 
again a PBE-like expression 11 j 



E, = j p{v) [er'(^.(r),C(r)) +J^(r.(r),C(r),<(r))] dr. 



with 



i7(r„C,t)-70'ln 1 + - 



7 1 + y4i2 + AH^ 



(5) 
(6) 



where p = p-^ -\- is the total density, Vs = [(47r/3)p]^/^ 
is the local Seitz radius, C, = (pf — Pi)/p is the relative 
spin polarization, = [(1 + C)^^^ + (1 ~ C)^^^]/2 is a spin 
scaling factor, e™'^ is the correlation energy per particle 
of the uniform electron gas, A is a function of e™'^ and 
(see Ref. |Tl|), t = (37r)i/^|Vp|/(4(;ip^/6) jg correlation 



reduced gradient, and 7 = (1 — In2)/7r2 is a parame- 
ter fixed by uniform scaling to the high-density limit of 
the (spin-unpolarized) correlation energy. The parame- 
ter f3, which is fixed to 0.066725 in the original PBE by 
the second-order gradient expansion of the correlation 
energy, is obtained for PBEint by fitting to a jellium sur- 
face reference system: /3 = /jPBEmt _ p. 052. In fact, the 
PBEint correlation can reproduce with high accuracy the 
wave-vector analysis of the correlation surface energies 
of jellium slabs of different thickness and rj [2^, outper- 
forming other PBE-like GGA functionals. This result is 
important to prove the ability of PBEint to describe ac- 
curately systems where different density regimes coexist. 
The value of /jP^Emt ^^so intermediate between those 
of PBE and PBEsol assuring an accurate compromise for 
a large palette of systems ranging from atoms to solids 

M. 
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Hybrid functional 

To construct our hybrid functional we follow the adia- 
batic connection scheme introduced in Ref. [sj where the 
XC (hybrid) functional is obtained as the result of the 
coupling-constant integration 

Elf = f Elt^{\)dX, (7) 

and the following ansatz is used for the hybrid-functional 
coupling-constant decomposition 

Elt'W = E^c^W + - - A)"-\ (8) 

with E^^ being the Hartree-Fock exchange energy, E^'^^ 
being the exchange energy of a given GGA functional, 
and E'^I^^{X) its coupling-constant decomposition. The 
parameter n controls the balance between the nonlocal 
Hartree-Fock and the local GGA exchange at different 
values of the coupling constant and is related to pertur- 
bation theory considerations [34 [ . 

If in Eq. ^ we use GGA=PBE and n=4, after per- 
forming the integration the PBEO functional 
(also known as PBEIPBE) is obtained. In this work 
instead we consider GGA^PBEint and we end up with 
the hPBEint fimctional 

1 



pihPBEint _ pPBEint 

XC XC 



(i?,"^-i?™). (9) 



In this case we do not fix the value of the parameter n, 
but we will consider n= 4 ,5, and 6 in order to assess 
the optimal value of this parameter on a broad range of 
situations. In fact, using a value of n = 4 can be better 
to describe molecular systems (as in PBEO), while larger 
values of the parameter might be more appropriate for 
transition-metals, clusters and interfaces. 



Computational details 

We tested the PBEint and hPBEint functionals for a 
series of molecular and solid state properties including 
atomization energies, structural properties, non-bonded 
interactions and reaction energies. The calculations were 
performed with the PBEint functional as well as with 
the hPBEint functional with n=4, 5, and 6. Moreover, 
calculations employing the PBE PBEsol and 
PBEO fsS*! functionals were also carried out for compar- 
ison. In all calculations a def2-TZVPP basis set H, 33| 
was employed. For transition metals the core electrons 



were replaced with effective core potentials (ECP) [38- 
l40j . All calculations were performed at the optimized 
ground-state geometry, except for those on the DM25 
and small interfaces test sets (see below). 

All calculations on molecular species were performed 
with a development version of the TURBOMOLE pro- 
gram package [4l'| . Calculations on bulk solids were per- 
formed with the FHI-AIMS program 42, 14^ using the 



light basis-set and a 18 x 18 x 18 fc-point grid. In this 
case, scalar relativistic effects were accounted for by the 
zeroth-order relativistic approximation (ZORA) [4J|. 

In more details, in the present work the following tests 
were considered: 

AE6: Atomization energies of Sill4 , SiO, S2, C3H4, 
C2II2O2, and C4H8; reference data were taken from 
Ref. 145|. Note that the AE6 test set was build 
to be representative for the results of the large 
Database/3 [i^, including 109 atomization ener- 
gies. 

W4: Atomization energies from the W4-08woMR 
test set of Ref. [i^l . It includes 83 atomization en- 
ergies of organic molecules selected from the origi- 
nal W4 test set [i^l excluding multi-reference cases. 

TMIO: Atomization energies of CrH, MnH, CoH, 
V2, Sc2, TiO, MnO, CuO, CrF, and CuF. Reference 
data were taken from Ref. 
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AUnAE: Atomization energies of the Au^, Au2, 
Au3, and Au5 clusters. Reference data, includ- 
ing relativistic and thermal corrections, were taken 
from Ref. [H 

K9: Barrier heights and reaction ener- 
gies of three organic reactions. Namely, 
0H-hCH4 ^CH3-hH20, H+0H->0-KH2, and 
H-f H2S— )-H2+HS. Reference data were obtained 
from Refs. S, \^ 

HB6: Binding energies of the hydrogen-bond inter- 
acting systems (H20)2, (HF)2, (NH3)2, NH3-H2O, 
(HCONH2)2, andlHC00H)2. Reference data were 
taken from Ref. |50|- 

DIG: Binding energies of the dipole-dipole inter- 
acting systems CH3CI-HCI, CH3SH-HCI, CH3SH- 
NCH, (H2S)2, (HC1)2, and H2S-HCI. Reference 
data were taken from Ref. (50|. 

MGBL19: Bond lengths of H2, CH4, NH3, H2O, 
HF, C2H2, HCN, H2CO, OH, CO, N2, F2, CO2, 
N2O, and CI2. Reference values were taken from 
Ref. 



51 



Note that this set provides a global as- 
sessment over bond lengths involving one hydro- 
gen atom and bond lengths not involving hydrogen 
atoms, which usually behave differently for differ- 
ent functionals (23| . 

AUnBL: Bond lengths of Auj", Au2, Au;^, Aug", 
Au4, Aug (capped pentagon, G^^ symmetry), 
SeAu2, and (ClAuPH3)2. Reference data were 
taken from Ref. [i^- 

F38: Harmonic vibrational frequencies of H2, CH4, 
NH3, H2O, HF, CO, N2, F2, C2H2, HCN, H2CO, 
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TABLE I: Mean absolute errors for different tests as obtained from PBE, PBEint, PBEsol, PBEO, and hPBEint calculations 
with n=4, 5, 6. For each line the best result is indicated by bold face, the worst is underlined . All energies are in kcal/mol 
(except for COH6). All bond lengths are in mA. Vibrational frequencies are in cm~^. 



Test 


PBE 


PBEint 


PBEsol PBEO 


n = 4 


hPBEint 

71 = 5 


71 = 6 








Organic molecules 










AE6 


14.5 


24.7 


34.9 


5.4 


12.4 


14.3 


15.7 


W4 


10.7 


15.5 


21.5 


5.7 


6.3 


7.5 


8.5 


OMRE 


6.8 


8.2 


12.0 


9.1 


12.8 


11.7 


11.1 


DC9 


10.6 


15.5 


17.6 


10.2 


13.6 


13.9 


14.2 


K9 


7.51 


9.09 


10.59 


3.96 


4.90 


5.62 


6.17 


MGBL19 


9 


10 


10 


6 


9 


8 


7 


F38 


56.8 


65.4 


65.9 


51.7 


53.2 


45.3 


45.6 








Transition metals 










TMIO 


13.4 


15.5 


18.3 


18.4 


16.2 


15.3 


12.0 


AUnAE" 


0.3 


2.2 


4A 


3.4 


2.5 


1.8 


1.3 


TMRE 


3.7 


6.9 


9.9 


11.1 


11.1 


9.4 


8.3 


AUnBL 


78 


31 


25 


86 


58 


56 


55 








Non-bonded interactions 








HB6 


0.4 


0.5 


M 


0.5 


0.7 


0.7 


0.6 


DI6 


0.4 


0.4 


LO 


0.4 


0.4 


0.4 


0.4 


LC6 
COH6'' 


59 
0.15 


32 
0.14 


Solid-state properties 
22 
0.21 









a) atomization energy per atom. 

b) eV/atom. 



CO2, N2O, CI2, and OH. Reference data were taken 
from Ref. 52 



OMRE: Reaction energies of organic molecules. It 
includes CH2+H2 ^ CH4, F2-I-H2 ^ 2HF, CeHg-h 
3H2 C6H12, CO + 3H2 ^ CH4 + H2O, SO2 + 
3F2 ^ SFe + O2, and C4H6 + C2H4 ^ CgHio. 
Reference data were taken from Ref. l53l 



TMRE: Reaction energies of transition-metal 
complexes. It includes Ni(C0)3 + CO Ni(C0)4, 
Fe(C0)4-hC0 ^ Fe(C0)5, iCl2 + CoCl2 ^ C0CI3, 
and 2FeCl2 — >■ Fe2Cl4. Reference data were taken 



from Ref. 21 



DC9: Nine reaction energies of medium-size 
molecules, which are usually treated poorly by DFT 
methods. Reference values were taken from Ref. 
47. 



LC6 and COH6: Lattice constants and cohesive 
energies of bulk Na (simple metal), Ag, Cu (transi- 
tion metals). Si, GaAs (semiconductors), and NaCl 



{ionic solid). Reference data were taken from Refs. 
il [HI The LC6 and C0H6 tests were not per- 
formed for hybrid functionals. 



DM25: Dipole moments of molecules. The test set 
considers a broad range of systems and of dipole 
moments, ranging from 0.11 to 11.56 Debye. In de- 
tail, it includes CO, CFCI3, furan, OCS, BF, CF2O, 
phenol, HCOOH, CH3SH, HNCO, CH3COOH, 
CH3OCHO, HF, H2O, CH3CI, CH3ONO, pyridine, 
H2C0, CICN, CuH, HON, CHOCH2OH, CH3NO2, 
LiCl, and NH2(CH=CH)6N02 (denoted as N6). 
Reference data and geometries of BF, H2O, CuH, 
H2CO, LiCl, and N6 are taken from Ref. [HB. The 
remaining reference data and (experimental) ge- 
ometries are taken from Ref. |57l All structures 
are available in supporting information [ssf . 



Small interfaces: Interaction energies of twelve 
small metal-molecule systems, representative of 
metal-molecule interfaces. The set includes Au2- 
SH, Au^-N2, AU2-SHCH3, Au3-SH, Au;]""N2, 
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TABLE II; Dipole moments (Debye) for different test molecules as computed with different DFT methods. The mean error 
(ME), mean absolute error (MAE) and mean absolute relative error (MARE) of each method with respect to reference values 
are reported in the last lines. For each row the best result is highlighted in bold style, the worst one is underlined. 



System 


PBE 


PBEint 


PBEsol 


PBEO 


n = 4 


hPBEint 

n = 5 


n = 6 


Rcf. 


CO 


0.22 


0.24 


0.23 


0.10 


0.11 


0.14 


0.15 


0.11 


CFCI3 


0.35 


0.33 


0.35 


0.42 


0.41 


0.40 


0.39 


0.45 


fur an 


0.59 


0.58 


0.59 


0.66 


0.65 


0.63 


0.62 


0.66 


OCS 


0.75 


0.74 


0.77 


0.80 


0.79 


0.78 


0.77 


0.71 


BF 


1.02 


1.04 


1.03 


0.94 


0.95 


0.97 


0.98 


0.79 


CF2O 


1.12 


1.13 


1.12 


1.12 


1.13 


1.13 


1.13 


0.95 


phenol 


1.25 


1.26 


1.26 


1.28 


1.29 


1.28 


1.28 


1.22 


HCOOH 


1.42 


1.42 


1.42 


1.47 


1.47 


1.46 


1.45 


1.41 


CH3SH 


1.58 


1.60 


1.61 


1.62 


1.64 


1.63 


1.62 


1.52 


HNCO 


1.99 


1.99 


2.01 


2.07 


2.07 


2.06 


2.05 


1.61 


CHqCOOH 

V.^ J. J. J V_y \_7 A.A. 


1.72 


1.72 


1.73 


1.79 


1.79 


1.77 


1.77 


1.70 


CH30CHO 


1.86 


1.86 


1.87 


1.89 


1.89 


1.88 


1.88 


1.77 


HF 


1.78 


1.79 


1.80 


1.83 


1.84 


1.83 


1.82 


1.82 


H2O 


1.87 


1.88 


1.89 


1.92 


1.93 


1.92 


1.91 


1.85 


CII3CI 


1.82 


1.83 


1.84 


1.89 


1.90 


1.89 


1.88 


1.87 


CII3ONO 


2.26 


2.27 


2.30 


2.21 


2.21 


2.23 


2.23 


2.05 


pyridine 


2.22 


2.23 


2.24 


2.27 


2.27 


2.26 


2.26 


2.19 


H2CO 


2.23 


2.24 


2.25 


2.42 


2.42 


2.38 


2.36 


2.39 


CICN 


2.96 


2.96 


2.98 


2.99 


2.99 


2.99 


2.98 


2.82 


CuH 


2.31 


2.33 


2.28 


2.85 


2.87 


2.76 


2.69 


2.97 


HCN 


2.93 


2.94 


2.95 


3.03 


3.04 


3.02 


3.01 


2.98 


CHOCH9OH 


2.40 


2.40 


2.43 


2.57 


2.57 


2.54 


2.51 


2.73 


CH3NO2 


3.41 


3.41 


3.43 


3.61 


3.61 


3.57 


3.54 


3.46 


LiCl 


6.90 


6.89 


6.90 


7.06 


7.06 


7.03 


7.01 


7.23 


N6 


16.94 


16.95 


17.03 


15.74 


15.77 


15.98 


16.14 


11.56 


ME 


0.20 


0.21 


0.22 


0.23 


0.23 


0.23 


0.22 




MAE 


0.35 


0.35 


0.36 


0.27 


0.27 


0.28 


0.29 




MARE 


13.5% 


14.4% 


14.4 


8.0% 


7.9% 


8.9% 


9.8% 




Std. Dev. 


1.10 


1.10 


1.11 


0.83 


0.84 


0.88 


0.92 





AU3-SCH3, Au4-SH, AuJ-N2, and AU4-SHCH3. 
The general structure of the molecule-cluster sys- 
tems was obtained from Refs. IsqI - IgiI . Each 
structure was then reoptimized at the TPSS/def2- 
TZVPP di, m [ei] level of theory. This geometry 
was successively employed in all the calculations. 
Reference values were obtained from CCSD(T) cal- 
culations 'gs'-'bs'] extrapolated to the complete basis 
set limit 66]. All structures are available in sup- 
porting information (Hsj . 

Model systems: We consider four sets of model 
systems for which exact or extremely accurate so- 
lutions are available. Namely, we consider: jellium 
surfaces, jellium clusters, Hooke's atoms, and the 



hydrogen atom with fractional spin. Additional de- 
tails are provided in the corresponding subsection 
later on. 



All tests are necessarily not exhaustive, because for 
computational reasons they only include a relatively 
small number of systems. Therefore, for specific cases 
they might be not fully representative of the true perfor- 
mance of the functionals. Nevertheless, since each test 
set was constructed to have good representativity and 
because the selected tests cover a broad range of proper- 
ties and systems, we believe that the present assessment 
offers a fair overview of the performance of the considered 
functionals. 
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TABLE III: Interaction energies (eV) of small organic molecules and gold clusters, computed with different DFT methods. 
The mean error (ME), mean absolute error (MAE) and mean absolute relative error (MARE) of each method with respect to 
CCSD(T) are reported in the last lines. For each row the result in best agreement with CCSD(T) calculations is highlighted 



in bold style, the 


one in worst 


agreement 


is underlined. 












System 


PBE 


PBEint 


PBEsol 


PBEO 




hPBEint 




CCSD(T) 












n = 4 


n — 5 


n = 6 




Au2-SH 


1.99 


2.20 


2.35 


1.57 


1.72 


1.81 


1.87 


1.98 


AU+-N2 


0.99 


1.13 


1.27 


0.84 


0.95 


0.99 


1.01 


1.30 


AU2-SHCH3 


0.88 


1.04 


1.17 


0.78 


0.90 


0.92 


0.94 


1.24 


Au3-SH 


3.02 


3.20 


3.33 


2.96 


3.09 


3.11 


3.12 


3.15 


AU3-N2 


0.77 


0.93 


1.06 


0.67 


0.78 


0.74 


0.77 


0.78 


AU3-SCH3 


2.73 


2.90 


3.03 


2.66 


2.79 


2.81 


2.82 


2.88 


Au4-SH 


2.18 


2.36 


2.49 


2.09 


2.23 


2.25 


2.27 


2.48 


Au^-N2 


0.67 


0.80 


0.92 


0.58 


0.68 


0.70 


0.72 


0.74 


AU4-SHCH3 


0.95 


1.11 


1.24 


0.93 


1.04 


1.06 


1.07 


1.05 


ME 


-0.16 


0.01 


0.14 


-0.28 


-0.16 


-0.13 


-0.11 




MAE 


0.16 


0.12 


0.16 


0.28 


0.16 


0.14 


0.12 




MARE 


10.61% 


8.84% 


13.08% 


18.89% 


10.30% 


9.16% 


7.72% 




Std Dev 


0.14 


0.14 


0.15 


0.15 


0.14 


0.13 


0.12 





RESULTS 

Quantum chemistry and solid-state benchmarks 

In Table U we report the mean absolute errors (MAE) 
for different tests, as obtained from PBEint and hPBEint 
calculations with n = 3, 4, 5, with the aim of assessing the 
performance of the functionals for a broad set of prob- 
lems. To this end, also PBE, PBEsol, and PBEO resuhs 
are reported, since these are natural references for the 
functionals considered here. Full results of all tests are 
available in supporting information [58| . 

For atomization energies of organic molecules, it is well 
known that GGAs are not very accurate (PBEsol largely 
fails), while the hybrid functionals perform better than 
the corresponding GGAs and this is more so for those 
functionals including a larger fraction of Hartree-Fock 
exchange Consequently, PBEO, with a MAE of 5 

kcal/mol, is the best functional overall and hPBEint with 
n = A yields the best performance among the different 
variants of hPBEint, reaching almost the PBEO accuracy 
for the W4 test. Note that even better results can be ob- 
tained by considering higher fractions of exact exchange, 
as shown for example in Ref. [2^ where the admixture of 
PBE GGA and 32% of exact exchange was found to yield 
the best results for thermochemistry. 

Whereas hybrid functionals improve significantly at- 
omization energies, for reaction energies involving or- 
ganic molecules the use of h ybr id functionals is not so 
systematically beneficial 47|, |67|, |68|. In fact, inclusion 
of Hartree-Fock exchange brings some improvements in 



some cases, e.g. for the K9 test where PBEO and hP- 
BEint with n = 4 perform best (see also the BH76RC 
test {29!, '47']), but has small effects for other cases, e.g. 
G2RC, NBPRC and DC9 tests Furthermore, in 

most situations improvements are smaller or even absent 
when large fractions of Hartree-Fock exchange are used; 
see for example BHLYP vs BLYP and B3LYP for the 
BH76RC test [li E^]. In particular, for the DC9 and 



OMRE benchmarks considered here, we find that PBEO 
brings almost no improvement or even a worsening of 
the results with respect to PBE. Analogously, the best 
results at the hPBEint level are obtained when n = 6 
and a relatively small fraction of nonlocal exchange is 
used. 

Different conclusions are found when one considers 
the energetic properties of transition-metal complexes or 
clusters (TMIO, AUnAE, TMRE). In this case in fact the 
use of PBEO corresponds always to a considerable wors- 
ening of the results with respect to PBE. A worsening of 
the performance is also observed in all cases for hPBEint 
with respect to PBEint when n = 4, although not so 
marked as for PBEO vs PBE. Improvements with respect 
to PBEint are instead found for hPBEint with n — 6, 
which becomes the best functional for TMIO, but not for 
the TMRE test. In this last case hPBEint with n = 6 is 
definitely the best hybrid functional in the present inves- 
tigation, but slightly worst than PBEint and more than 
twice worst than PBE. Note however that the TMRE test 
is a particular difficult test in the context of the present 
work, because it requires a delicate balancing between 
the description of metal complexes (best described by 
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PBEint and PBE) and small organic molecules (best de- 
scribed by hybrids with a large fraction of Hartree-Fock 
exchange) . 

Concerning structural properties, hybrid functionals 
provide a modest improvement for organic molecules 
(MGBL19 and F38), while they worsen the results for 
gold clusters (AUnBL). Interestingly, the hPBEint func- 
tional with n = 6 yields the best results among the differ- 
ent variants of hPBEint, outperforms PBEO for AUnBL 
and F38, and gives almost the same result as PBEO for 
MGBL19. Thus, the combination of PBEint with 1/6 of 
Hartree-Fock exchange proves to be a very good choice 
for structural properties of finite systems. 

For non-bonded interactions very good results are ob- 
tained at the GGA level using the PBE and PBEint 
functionals, while poorer results (especially for hydro- 
gen bonds) are given by the PBEsol functional. Hybrid 
functionals cannot improve the performance of PBE and 
PBEint, yielding the same MAE as the GGAs for the DI6 
test and slightly larger errors for the HB6 test. In this 
latter case, the best results at the hybrid level are ob- 
tained by PBEO, taking advantage of the high accuracy 
of PBE for the hydrogen bonds, and by hPBEint with 
n=6, taking advantage of the small fraction of Hartree- 
Fock exchange included. 

Concerning solid-state properties, PBEsol is the most 
accurate for lattice constants, while PBEint is interme- 
diated between it and PBE. However, the limited accu- 
racy of PBEsol for atomic energies causes it to perform 
poorly for cohesive energies (C0H6 test). In this case, 
where a good description of both bulk solids and atoms is 
required, PBEint shows a remarkably good performance, 
due to its well balanced description of different density 
regimes. 

Dipole moments 

In Tab. |TT] we report the performance of the different 
functionals for the description of dipole moments of sev- 
eral test molecules. Benchmarking the dipole moments 
is not as common as energy and structural tests, never- 
theless it is very important to verify the accuracy in the 
description of density [g^ [t^ . This test allows in fact to 
probe indirectly the quality of the description, provided 
by different functionals, of the electron density distribu- 
tion in molecules, especially in valence and asymptotic 
regions. This issue is particularly relevant for GGA func- 
tionals, because due to the Coulomb self-interaction er- 
ror, often they yield a poor description of the electron 
density far from the atomic core. For the same reasons, 
an improved description is usually obtained by hybrid 
approaches. 

The general trends outlined above are confirmed by 
inspection of the results of Tab. [Hi The smallest er- 
rors are obtained by the hybrid functionals, which yield 



a global MAE of about 0.28 Debye. A MAE larger by 
30% is found instead for GGA methods which perform all 
very similarly and display a global mean absolute error of 
about 0.35 Debye. In particular, the best results are ob- 
tained when the larger fraction of Hartree-Fock exchange 
(0.25 for PBEO and hPBEint with n=4) is considered, 
having a MAE of 0.27 Debye and a MARE below 8%. 
Nevertheless, good results are achieved also by hPBEint 
with n=6, including only a relatively small fraction of 
Hartree-Fock exchange. 

We note however that the results of Tab. |ll] show for 
all the functionals several important deviations from the 
trends depicted above. In fact, for some systems (e.g. 
HNCO, pyridine, CICN) the best results are achieved at 
the GGA level, while PBEO yields the worst agreement 
with the reference. Thus, none of the functionals con- 
sidered here can be labeled as fully reliable for the cal- 
culation of the dipole moment of an arbitrary molecular 
system, but the best compromises are provided by the 
hPBEint functional with n=6, which yields the largest 
number of results with errors <0.1 Debye, and PBEO 
that yields the best overall MAE and MARE and the 
largest number of results with errors <0.2 Debye. 

Metal-molecule interfaces 

In this subsection we investigate the ability of different 
functionals to describe hybrid metal-molecule interfaces. 
These systems pose a difficult challenge to any compu- 
tational method because of the very different theoreti- 
cal and numerical issues raised by metallic (extended) 
systems and molecules. The former are in fact mainly 
characterized by a slowly-varying density regime and 
can be properly described with local or semilocal ap- 
proaches. The latter instead display a significant contri- 
bution from rapidly-varying density regions and require 
a proper treatment of nonlocal interactions. 

In view of practical applications it would be nice, of 
course, to perform tests for the interaction of (function- 
alized) organic molecules with relatively large metallic 
clusters or extended (eventually semi-infinite) surfaces. 
However, for these systems no benchmark value to assess 
the performance of the methods can be obtained. Accu- 
rate experimental values are in fact extremely rare, while 
high level approaches cannot be applied due to their ex- 
ceedingly high computational cost. Note also that for 
large metallic systems static and high-order correlation 
effects start to play a fundamental role, thus accurate 
results can only be obtained through very sophisticated 
methods. 

For the reasons mentioned above, in this paper we re- 
strict our attention to a set of relatively small metal- 
molecule systems composed of small organic molecules 
interacting with gold clusters of 2, 3 and 4 atoms. For 
such systems quite accurate benchmark results have been 
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TABLE IV: Mean absolute errors (MAE) in Hartree or mean absolute relative errors (MARE) for different functionals for the 
proposed model systems. Best results are in bold style; worst results are underlined. 



System 


PBE 


PBEint 


PBEsol 


PBEO 




hPBEint 














n = 4 


n — 5 


n = 6 


Jellium surf, energy (MARE) 


3.15 


2.29 


2.32 


1.88 


2.29 


2.29 


2.29 


Jellium spheres (MAE) 


0.46 


0.46 


0.56 


0.57 


0.23 


0.18 


0.17 


Hooke's atom (MARE) 


3.63 


3.45 


3.41 


3.87 


3.78 


3.71 


3.66 


Fractional spin H (MAE) 


29.59 


27.78 


27.63 


42.31 


40.80 


38.19 


36.46 



TABLE V: XC jellium surface energies (erg/cm^) of the global hybrids PBEO and hPBEint. The fixed-node diffusion Monte 
Carlo (DMC) calculations ^731, III] are reported for reference. (Ihartree/bohr^ = 1.557 x lO^erg/cm^.). PBE, PBEint, and 
PBEsol results are taken form Table I of Ref. 





PBE 


PBEint 


PBEsol 


PBEO 


n = 4 


hPBEint 
n — 5 


71 = 6 


DMC 


2 


3265 


3378 


3374 


3312 


3378 


3378 


3378 


3392 ± 50 


3 


741 


774 


774 


756 


774 


774 


774 


768 ± 10 


4 


252 


267 


267 


259 


267 


267 


267 


261 ± 8 


6 


52 


56 


56 


55 


56 


56 


56 


53 ±... 



obtained with coupled cluster single, doubles with per- 
turbative triple correction (CCSD(T)) calculations and 
an assessment of DFT methods could be carried out. 

Results are reported in Tab. IIIIl At the GGA level, 
PBE presents a general tendency to underestimate the 
interaction energies, while a slight overestimation is ob- 
served by PBEint. Finally, PBEsol shows a rather 
strong overestimation of the interaction energies. Over- 
all, PBEint (with a MAE of 0.16 kcal/mol) is slightly 
better than PBE and PBEsol. In addition, we note that 
most of the accuracy for PBE comes from the smallest 
systems (e.g., Au2-SH) and larger errors are found when 
the system size increases. For PBEint an opposite be- 
havior is observed. This trend was already observed for 
gold clusters ^22] and led to the conclusion that PBEint 
provides the most accurate description of medium and 
large metal clusters. A similar extrapolation can be sup- 
posed to apply also in the present case of metal-molecule 
interfaces, suggesting that superior results can probably 
be achieved by PBEint for realistic interfaces. This con- 
clusion is also supported by the observation that for in- 
terfaces involving large systems slowly-varying density 
regions will gain a more important role. Thus PBEint 
(and partially PBEsol) will be favored over PBE. 

The addition of a fraction of Hartree-Fock exchange in 
the functionals has the effect to reduce the computed in- 
teraction energies. This results in a worsening of PBEO 
with respect to PBE, but an improvement of hPBEint 
with respect to PBEint. In fact, all the variants of hP- 
BEint show a MARE lower than PBE and PBEint. In 
particular, the option with n = 6 displays the smallest 



MAE and MARE of all the functionals considered in this 
work. We note also that, because of the small fraction 
of Hartree-Fock exchange that it includes, hPBEint with 
n = 6 may be supposed to be the best compromise for 
interfaces of interest in practical application, where large 
metallic systems are involved 7l|, |72| . 



Model systems 

In this subsection, we show the performance of PBEO 
and hPBEint for several important model-systems: 

(i) jellium surfaces, that contain the main physics of 
simple metal real surfaces; 

(ii) jellium spheres, that are simple models for simple- 
metal-particles; 

(iii) Hooke's atom, that represents two interacting elec- 
trons in an isotropic harmonic potential of frequency oj. 
At small values of w, the electrons are strongly correlated, 
and at large values of oj, they are tightly bounded, two 
important cases in many condensed matter applications. 

(iv) The hydrogen atom with fractional spin, a model 
for the static correlation given by degenerate states. 

For all these cases, we use non-self-consistent accurate 
calculations (for jellium surfaces we use numerical LDA 
orbitals and densities, for jellium spheres we use accurate 
exact exchange orbitals and densities, for the Hooke's 
atom we use exact orbitals and densities, and for the spin- 
dependent hydrogen atom we use exact-exchange orbitals 
and densities.) For the hybrids, we use exact-exchange 
instead of the Hartree-Fock one. 
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The mean absolute (relative) errors of different func- 
tionals for the selected model problems are summarized 
in Tab. |TVl A detailed analysis of the results is given in 
the following. 



Jellium surfaces 

The PBEint GGA performs remarkably well for the 
jellium surfaces, giving very accurate exchange and XC 



0.001 



surface energies [24 



26j . Because there are no error com- 
pensations between the exchange and correlation parts 
and the PBEint exchange surface energies are extremely 
close to the exact-exchange ones j23| , any PBEint hybrid 
will yield the same results as the PBEint GGA and thus 
be very accurate for the jellium surfaces, as we explicitly 
show in the Table fVl This is a very important result for a 
functional designed for hybrid interfaces, where the sur- 
face physics plays a dominant role. We also mention that 
PBEO improves considerably over PBE, but still it is not 
so accurate in the range 2 < < 3, where important 
metals lie (e.g. Al has Vg = 2.07, and Cu has Ts = 2.67). 

We mention that the RGE2 GGA of Ref. [Z5| gives 
similar jellium xc and x-only surface energies as PBEint 
(see Table 3 and Fig. 3 of Ref. [75] for a detailed analysis 
of jellium surfaces). Both RGE2 and PBEint have the ex- 
change enhancement factor Fx — > 1 + ^s^ + O^s^) at small 
s, and thus the requirement of vanishing O(s^) terms in 
this expansion seems very important for accurate jellium 
x-only surface energies at the GGA level. Note that the 
exchange part of revTPSS meta-GGA [t^, which recov- 
ers the 4-th order gradient expansion of the exchange 
energy, is also accurate for jellium surfaces. 



Jellium spheres 

In Fig. [l] we show the errors in total energies per 
electron {E'''='"'yN - E'^pp''"'-' /N) versus the number of 
electrons N, for jellium spheres of magic numbers (iV=2, 
8, 18, 20, 34, 40, 58, 92, and 106) with bulk parameter 
rg = 4. (Note that sodium has = 3.93.) The refer- 
ence values E'^^"-'^^ /N represent the corrected DMC data 



(See Eq. (40) and Fig. 5 of Ref. iTTl). We note that 
HF and EXX perform similarly and badly for this prob- 
lem [tI], whereas LDA is remarkably accurate for large 
clusters [t^ [s^. We also recall that self-consistent ef- 
fects are very small (see Table I of (sij), such that our 
non-self-consistent results are expected to give an hon- 
est and accurate picture of the performance of various 
functionals for jellium (metal) particles. 

We observe that PBEO does not improve upon PBE, 
being even worse for N < 34. On the other hand, hP- 
BEint (with n=4, 5 and 6) significantly improve over the 
original PBEint for N < 20. However, on average the 
best performance is given by hPBEint (n=6). This is 



-0.001 



-0.002 
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FIG. 1: Errors in total energies per electron {E'^^"-'^*' jN — 
£,approx^^. Hartree) of jellium spheres for = 4. E"""""^ jN 
represents the corrected DMC data [t^I • We use magic neutral 
clusters with N=2e~ , 8e", 18e", 20e", 34e", 40e", 58e", 
92e", and 106e". 



also a very important result for a functional designed for 
hybrid interfaces, where metal particles play a crucial 
role. 



Hooke's atom 

In Fig. [21 we report the relative absolute errors of the 
hybrids, for the Hooke's atom with several frequencies w. 
In the tightly-bounded regime {tq = (a;^/2)~^/'^ is small) 
all the functionals performs similarly, though PBEO is 
the best in this extreme case. On the other hand, for 
To > 5, that is the physical case for most of the real ap- 
plications, the hPBEint hybrids are always better than 
PBEO, with a superior performance for the n = 6 case. 
However, the errors increase when the strong correlation 
increases, showing that these global hybrids can not de- 
scribe satisfactorily the strongly correlated regime. 



Hydrogen atom with fractional spin 

For the static correlation case, the exact exchange fails 
badly, and thus the global hybrids usually increase the 
error in the static correlation with respect to the orig- 
inal GGA. So, from this point of view, the HF mixing 
should be as small as possible. In Fig. [31 we show that 
indeed for the different variants of hPBEint increasingly 
good results are obtained in the series n=4, 5, 6. Thus, 
hPBEint with n = 6 has the best performance, while the 
worst results are obtained with PBEO, which is slightly 
worst than hPBEint with n=4. 

We recall however that the iV-electron self-interaction 
error 82|, or the delocalization error [s^l, that is related 



to the convexity behavior of the functional at fractional 
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0.03 



0.02, 



0.01 




FIG. 2: AS,c = KE^PP - (Hartree) for the 

Hooke's atom with difTerent frequencies, ro = {uj'^ /2)~^^^ is 
the classical electron distance 12611. 
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FIG. 3: AE = E[n,C] ~ E[n,( = 1] (Hartree) versus the 
spin-polarization (. 



number of electrons, decreases when the mixing parame- 
ter of a global hybrid increases. Thus, the errors of static 
correlation and of delocalization can not be both reduced 
by a global hybrid. 



suits, the benchmarks were divided in two subsets: the 
molecular-based tests including organic molecules, non- 
bonded interactions, transition metals and dipoles and 
the other-systems including interfaces, solid-state and 
model systems. For the two groups RMAE-mol and 
RMAE-other indicate respectively the average RMAE. 
Finally, RMAE-tot indicates the averaged RMAE among 
all seven benchmarks. 

Among the considered GGAs, PBEint shows a remark- 
ably good performance yielding the best global RMAE- 
other and RMAE-tot and a RMAE-mol not much higher 
than PBE and even lower than PBEO. In fact, despite 
PBEint provides the smallest MAE only in few cases, it 
has a more balanced performance on different tests. On 
the contrary, PBEO and PBE give very accurate results 
in some cases but behave rather poorly for other tests. 
The good balance over a broad set of systems and prop- 
erties is an important feature of the PBEint functional 
which makes it suitable for applications in complex sys- 
tems where different situations (or density regimes) may 
coexist. 

Concerning the hybrid functionals similar considera- 
tions apply. We observe that functionals incorporating 
25% of Hartree- Fock exchange (PBEO and hPBEint with 
n = 4), despite working very well for energetic proper- 
ties of organic molecules (see Table H]) , perform globally 
worst than their GGA counterparts and other hybrids 
with a lower Hartree-Fock content, again because of im- 
portant failures for some properties (e.g., structures and 
reaction energies). The best overall results are obtained 
with hPBEint using n — Q, which gives also a RMAE- 
mol of 0.95, (better than PBEint and PBEO). Notably, 
hPBEint n = 6 has in fact a nicely uniform performance 
for all the tests considered in the present work, being 
never the worst one and showing a good accuracy for all 
systems and properties. 

Moreover, hPBEint with n = 6 and PBEint are the 
only functionals that can describe rather accurately at 
the same time energetic and structural properties (this 
holds also for bulk in the case of PBEint). 



CONCLUSIONS 



Global assessment 

To compare the MAE of different tests, we computed, 
for each group of tests, the MAE relative to PBEint de- 
fined as 



RMAE (method) 



M AEi (method) 
^ MAi;, (PBEint)' ^ ' 



where i runs over a given set of tests. The RMAEs for 
the different classes of tests considered in this work are 
collected in Tab. I VII To help the analysis of the re- 



in this paper, we performed an extended study of the 
PBEint global hybrids. We tested important energeti- 
cal and structural properties of organic molecules and 
metal complexes, as well as metal-molecule interfaces. 
We showed that the hPBEint functional, with the opti- 
mal mixing parameter n=6, maintains the well-balanced 
behavior of the PBEint GGA functional, improving the 
overall accuracy, over a wide range of problems and can 
thus be successfully employed in complex studies such as 
the description of hybrid metal-molecule interfaces. 

In fact, the PBEint functional and its hybrid forms 
take benefit from the good treatment of both the slowly- 
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TABLE VI: Mean absolute errors relative to PBEint (RMAEs; see Eq. (|10p ) for different sets of tests and all functionals. 
In the last lines are reported: the global RMAE-mol for molecular properties (organic molecules, non-bonded interactions, 
transition-metal systems, dipoles); the global RMAE-other for properties not concerning molecular systems (interfaces, solid- 
state systems, model systems), in this case for hybrid functionals solid-state RMAEs are not considered; the global RMAE-tot 
for all properties. For each row the best results are highlighted with bold style, the worst ones are underlined 



Set 


PEE 


PBEint 


PBEsol 


PBEO 


n = 4 


hPBEint 
n = 5 


n = 6 


Organic molecules 


0.76 


1.00 


1.21 


0.60 


0.79 


0.78 


0.79 


Non-bonded inter. 


0.81 


1.00 


2.84 


0.93 


1.15 


1.11 


1.09 


Transition metals 


1.01 


1.00 


1.35 


1.84 


1.41 


1.24 


1.09 


Dipoles 


0.99 


1.00 


1.01 


0.76 


0.77 


0.80 


0.82 


RMAE-mol 


0.85 


1.00 


1.50 


1.03 


1.04 


0.97 


0.93 


Interfaces 


1.09 


1.00 


1.19 


1.30 


1.03 


0.98 


0.93 


Solid-state 


1.46 


1.00 


1.09 










Model systems 


1.12 


1.00 


1.05 


1.18 


1.02 


0.96 


0.93 



RMAE-other 


1.22 


1.00 


1.11 


1.24° 


1.03° 


0.97° 


0.93° 


RMAE-tot 


1.03 


1.00 


1.39 


1.10° 


1.03° 


0.98° 


0.94° 



a) It does not include solid-state. 



and rapidly- varying density regimes, granted by Eq. 
Therefore, despite some limitations for the thermochem- 
istry of organic molecules, they provide reasonably accu- 
rate results for all the tests considered and outperform 
the other functionals for transition metals, hybrid in- 
terfaces and solid-state systems. For these reasons the 
PBEint and hPBEint(n=6) functional are especially ap- 
pealing for applications involving transition-metal clus- 
ters and/or the interaction of metal clusters with organic 
molecules. On the other hand, because Eq. (j4]) is based 
on a simple GGA ansatz, PBEint-based functionals can- 
not be expected to outperform specialized functionals for 
specific problems (e.g. atomization energies of small or- 
ganic molecules). 

In the case of hPBEint the good behavior of the under- 
lying PBEint GGA allows to achieve a good compromise 
in the overall performance using only a small fraction of 
Hartree-Fock exchange (16.67% for n—6), similarly with 
the case of hybrid meta-GGAs. This suggests that the 
PBEint GGA functional could be an interesting start- 
ing point for the construction of more advanced hybrids 
(e.g. local hybrids or screened ones) of good accuracy 
and large applicability in computational studies. 

Finally, we note that a recent work [s^l showed that a 
GGA can be accurate for both solids and molecules, if it 
satisfies a statistical constraint for atomization energies 
derived from simple one-electron densitites. In particu- 
lar, the zPBEsol and zPBEint functionals of Ref. [13] 
significantly improve the atomization energies and other 
spin-dependent properties, performing the same as the 
original functionals for any closed-shell system. Thus, 



they can be a good starting point for future developments 
of the hPBEint functional. 
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